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We study the mechanics and statistical physics of dislocations interacting on cylinders, motivated 
by the elongation of rod-shaped bacterial cell walls and cylindrical assemblies of colloidal particles 
subject to external stresses. The interaction energy and forces between dislocations are solved 
analytically, and analyzed asymptotically. The results of continuum elastic theory agree well with 
numerical simulations on finite lattices even for relatively small systems. Isolated dislocations on 
a cylinder act like grain boundaries. With colloidal crystals in mind, we show that saddle points 
are created by a Peach-Koehler force on the dislocations in the circumferential direction, causing 
dislocation pairs to unbind. The thermal nucleation rate of dislocation unbinding is calculated, for 
an arbitrary mobility tensor and external stress, including the case of a twist-induced Peach-Koehler 
force along the cylinder axis. Surprisingly rich phenomena arise for dislocations on cylinders, despite 
their vanishing Gaussian curvature. 
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INTRODUCTION 

Defects in crystals such as dislocations have been stud- 
ied extensively for more than seven decades [IH2], and 
their importance in condensed matter physics and ma- 
terial science is widely recognized jl]. Systems of dislo- 
cations in both two and three dimensions can be real- 
ized experimentally. Mechanical properties of bulk met- 
als are strongly affected by the dynamics of dislocation 
lines within them [4], and two-stage melting of a two- 
dimensional crystal can be driven by dislocations [5] . In- 
teresting applications involve a two dimensional particle 
array with a periodic boundary conditions in one direc- 
tion: a cylindrical crystal. One such example involves 
interacting colloids on the surface of a liquid film coat- 
ing a solid cylinder, where repulsive forces give rise to 
the self-organized emergence of a two dimensional crys- 
talline solid. Defects in colloidal assemblies on the related 
curved surfaces of capillary bridges were recently studied 
experimentally and theoretically [6]. Here, the Gaussian 
curvature can be positive or negative; the zero Gaussian 
curvature of a cylinder is a special case [7J. The growth 
of the cell walls of rod-shaped bacteria provides a bio- 
physical example. Their geometry can be approximately 
described by a cylinder, and in a recent study [8l |9] we 
have argued that one may regard cell wall elongation as 
mediated by dislocation climb [1]. As a final motivation 
we note that the hydrodynamic interactions of vortices 
on cylinders is mathematically similar to those of dislo- 
cations, albeit simpler since in this case the bare inter- 
actions on a flat surface are isotropic [10] . Vortices on 
superfluid Helium films were studied in a cylindrical ge- 
ometry, in order to model superfluids in porous materials 
[11) . A type II superconducting, hollow wire could have 
similar interacting vortices on its cylindrical surface. 

Fig. ([I]) illustrates the simplest dislocation in two- 
dimensional square and triangular lattices. For the 




FIG. 1. Lattice configuration corresponding to a single dis- 
location with Burgers vector in the positive x direction for 
a square (left) and triangular (right) lattice of masses and 
harmonic springs relaxed to their equilibrium configuration. 
The Burgers vector is defined as the deficit (red arrow) of a 
clockwise loop that would be closed in the absence of a dis- 
location; hence b = bx in both cases, where b is the lattice 
constant. The 'T' symbol encodes both the position and ori- 
entation of the dislocation; the leg of the inverted 'T' points 
in the direction of added material. 



square lattice, the dislocation can be thought of as the 
termination of a semi-infinite row of lattice points added 
to the upper half of an otherwise perfect lattice. For the 
triangular lattice, the dislocation represents the termi- 
nation of two semi-infinite rows of particles, inclined at 
±30° to the vertical. Generally, the dislocation can be 
characterized by a topological invariant called the Burg- 
ers vector, defined by the Burgers circuit [I], b = § jjfdl 
(see Fig. [I]). It is immediately apparent that such de- 
fects distort the crystalline lattice in their vicinity. For 
both the triangular and square lattice, it is possible to 
associated a particular set of neighbors to each lattice 
site by assigning near neighbor bonds. This allows us 
to define the dislocations using the number of neighbors: 
in the square lattice, the dislocation core has three sites 
with five nearest neighbors. In the triangular lattice, 
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the Voronoi construction shows that the dislocation is 
characterized by a five-fold site (a site with five nearest 
neighbors) bonded to a seven-fold site. 

The structure of the manuscript is as follows. We first 
describe two scenarios, one motivated by cell wall growth 
of rod-shaped bacteria and the other by colloidal crystals 
on cylinders, and show that they require understanding 
the dynamics of interacting dislocations on a rectangu- 
lar and triangular lattice, respectively. In contrast to 
the predominantly climb dislocation dynamics relevant 
to elongating bacteria [9], the dynamics appropriate to 
colloidal assemblies on cylinders is predominantly glide, 
i.e., the motion is parallel as opposed to perpendicular 
to the Burgers vector [3]. We analyze the form of the 
interactions between dislocations with arbitrary Burgers 
vector and compare the results to numerical simulations 
showing that, surprisingly, good agreement with contin- 
uum elastic theory is already achieved for relatively small 
systems. We then proceed to exploit a useful connec- 
tion between isolated dislocations on a cylinder and the 
physics of grain boundaries. Finally, we discuss related 
'nucleation' problems: Upon the addition of a force in 
the circumferential direction driving defects of opposite 
sign apart, there will be a finite unbinding rate at non- 
zero temperature. We calculate this rate using Langer's 
generalization of Kramers' theory [T2 H2], and find inter- 
esting geometrical effects associated with the cylindrical 
geometry. We also discuss the effect of a twisting stress 
applied to the ends of a cylinder coated with colloids; in 
this case the strain is relaxed by dislocation pairs sepa- 
rating predominantly along the cylinder axis. The Airy 
stress function for a dislocation on a cylinder is calcu- 
lated in Appendix A. In Appendix B, we discuss some 
subtle aspects of the quantization of stresses and strains 
on a cylinder, due to periodic boundary conditions. 



Bacterial cell wall growth 

Bacterial cell walls are made of a partly ordered mesh 
of peptidoglycan [TJ] , which can be only a single molecule 
layer thick in gram-negative bacteria. While there are 
still many open questions regarding architecture and 
growth, this meshwork is known to consist of circumfer- 
ential glycan strands cross-linked by peptides, as shown 
schematically in Fig. [2j We note that rod-shaped bacte- 
ria with very large aspect ratios can be created by sup- 
pressing the septation process associated with cell divi- 
sion [15] . To insert new material into the structure, de- 
fects in the mesh have to be created. Inserting a single 
glycan strand between two existing ones would not pre- 
serve the topology of the network. However, inserting two 
glycan strands between two existing ones does preserve 
the topology. This observation has led to the "three-for- 
one" hypothesis, see for example Ref [IB], or, alterna- 
tively, the more symmetric process shown in Fig. [2j Each 
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FIG. 2. Schematic illustration of insertions of new glycan 
strands into the peptidoglycan mesh of a bacterial cell wall. 
The new material is inserted in the vicinity of the blue ar- 
row, which points azimuthally around the cylinder and sits 
near the core of the dislocation. The axis of the cylinder runs 
horizontally, which is also the direction of the Burgers vec- 
tor of this dislocation. Although the the region around the 
new insertions is distorted, the connectivity of the structure 
is preserved locally due to the simultaneous insertions of two 
extra glycan strands into the structure, respecting the lattice 
geometry. The rectangular (dashed-line) box shows the bio- 
logically relevant unit cell, which is however not the minimal 
unit cell of the underlying lattice [Sj. 



zig-zag glycan strand is composed of alternating sugar 
units called NAM (N-acetylmuramic acid) and NAG (N- 
acetylglucosamine) . For the unit cell shown in Fig. [2] 
inserting two glycan strands is equivalent to the addi- 
tion of a single additional unit cell. The coarse-grained 
lattice obtained in this way is rectangular, and the dislo- 
cation can be mediated by edge dislocations climbing in 
this lattice in the circumferential direction. The relevant 
Burgers vectors point in the direction of the cylinders 
axis of symmetry [8l [9] . 

Here, we shall simplify the analysis by neglecting the 
anisotropy associated with the two-dimensional rectan- 
gular lattice, which would require a non-isotropic elastic- 
ity theory with four elastic coefficients to describe elas- 
tic deformations of the structure [17]. Rather, we shall 
approximate the systems's free energy by the standard 
isotropic form |18] : 



F, 



elastic 



(1) 



where Uij(x) — \ [diUj(x) + djUi(x)] is the 2d strain ten- 
sor. A and fi are the two-dimensional Lam coefficients. 
Fig. |3]shows a number of dislocations in a square lattice. 
Dislocation climb, mediated by glycan strand extension 
machinery, can be affected by interactions between these 
dislocations. 
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FIG. 3. Schematic of dislocations (i.e. glycan strand ends) 
on the cylindrical portion of a bacterial cell wall of radius Rq. 
A subset of these defects rotate circumferentially when pro- 
pelled by the addition of material from inside the bacterium, 
mediated by strand elongation machinery on the dislocation 
cores (not shown). Constant velocity climb motion of such 
dislocations leads to exponential elongation of the cylinder 
length L(t); see Refs. 015]. 

Colloidal crystals 

Recent experimental advances allow the creation of 
dislocations in colloidal assemblies with fascinating in- 
terfacial geometries, allowing study of the interplay of 
geometry, including Gaussian curvature, with defects in 
the lattice structure [6] [19] . In this case the lattice of col- 
loidal particles is typically triangular, with lattice vectors 
that interact only weakly with the directions of principal 
curvature. Upon approximating the pair interaction be- 
tween colloids as harmonic with spring constant fc s for 
small displacements about the equilibrium positions, the 
lattice can be described elastically with isotropic effective 
Lam coefficients A = /i = \/%k s /& [20 . Fig. (HJ) shows an 
example of a single dislocation in a triangular lattice. As 
we shall show here, even in the absence of Gaussian cur- 
vature, the periodic boundary conditions associated with 
a cylindrical geometry give rise to novel phenomena, not 
found in an infinite plane (to be referred to as "flat space" 
in the following). In flat space, rotational invariance al- 
lows an arbitrary orientation of the crystallographic axis. 
On a cylinder, however, square and triangular lattices 
can have an energetically preferred orientation relative 
to the long axis of the cylinder [3T]. With bacterial cell 
walls, for example, it may be easier to bend the glycan 
strands than the alternating amino acid cross-bridges, 
which would lead to a preferred orientation on a cylin- 
der. In this paper, we shall focus primarily on triangular 
lattices with Bragg rows that run either azimuthally or 
along the cylinder axis. 

Phyllotaxis 

In various biological as well as non-living systems, sub- 
units are arranged in an ordered lattice, wrapped on 
a cylinder, a particular case of a phenomenon known 
as phyllotaxis (meaning "leaf arrangement" in ancient 
Greek). These range from plants [23] [53], rod-shaped 




FIG. 4. Model of a dislocation pair on a cylinder. This struc- 
ture can systematically elongate if new particles (green) are 
added by two counter-rotating dislocations moving via climb. 
For one of the dislocations, the climb direction is indicated by 
a blue arrow. For the other, the Burgers circuit around the 
dislocation is shown (yellow) as well as the resulting Burgers 
vector b. Alternatively, the dislocation pair could separate 
by glide (motion parallel to b) , with motion predominately 
along the cylinder, with negligible elongation. In both cases, 
however, a triangular lattice with slightly tilted Bragg rows is 
created in the center, relative to the purely azimuthal Bragg 
rows made of silver particles on the two sides. Although the 
model was constructed with magnetic beads interacting via 
dipole-dipole interactions, we expect similar configurations for 
micron-sized colloids assembled on a cylindrical substrate. 

viruses and bacterial flagella [24-26 to systems where 
subunits arrange via magnetic interactions [2TH2"9"] . Each 
phyllotactic arrangement can be characterized by two in- 
tegers (M,N), such that Mbe x + Nbe 2 = Wy, with e x 
and ei the two lattice vectors of the triangular lattice, 
forming a 60° angle between them, and y points in the 
azimuthal direction around the cylinder of circumference 
W = 2ttR. Such tessellations where recently found useful 
also in the context of the Thompson problem on a cylin- 
der, i.e., how colloids pack in the bulk of a cylinder |30) . 
In various scenarios the energetically preferred tessella- 
tion can depend on certain external conditions [31] [32] , 
and it is interesting to understand the dynamics of the 
process through which this change comes about. It is 
plausible that the boundary region between two compet- 
ing phyllotactic tessellations will consist of one or several 
dislocations, see for example Fig [4] as we shall show in 
section [] a dislocation on a cylinder is equivalent to a 
grain boundary. As an example, Ref [5S] suggests that 
such a change in the phyllotactic arrangement of the tail- 
sheath of the bacteriophage T4 is driven by 6 dislocations 
symmetrically arranged on the circumference of the cylin- 
der. These observations provide additional motivation 
for understanding the interactions between dislocations 
on a cylinder, which is a necessary step to quantify the 
"dynamical phyllotaxis" problem sketched above. 



INTERACTING DISLOCATIONS 

Consider an edge dislocation on the surface of a cylin- 
der, with radius R and infinite length. We denote the 
coordinate along the symmetry axis by x, and the other 
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by y, so that y is periodic with a period of W = 27ri?. 
For simplicity, and without loss of generality, let us de- 
fine the axes' origin at the dislocation core. In Ref. [9], 
some of the components of the stress tensor produced 
by edge dislocation with Burgers vector b = ±bx were 
evaluated. Here, we shall calculate all components of the 
stress produced by dislocations with a Burgers vector b 
either in the ±bx or ±by direction. For a lattice of any 
symmetry or orientation with respect to the cylinder axis 
we could always decompose the dislocation's Burgers vec- 
tor into two orthogonal components in these directions. 
Within linear elasticity (and assuming isotropic elastic 
constants) one could thus use the results presented be- 
low to find the stress field induced by a dislocation with 
a Burgers vector of arbitrary orientation. 

In flat space, one may consider a single dislocation in 
an otherwise perfect crystal, in which case the energy 
of the system diverges logarithmically with the system 
size. However, the energy of a single dislocation on a 
long cylinder of length L, in general diverges as the area 
WL, as the subsequent discussion of long-range strains 
will make clear. An important special case is a Burgers 
vector in the x direction, when the divergence is only 
logarithmic in W and independent of L. In the case of 
multiple dislocations on a cylinder, the energy can be 
decomposed into the divergent terms mentioned above, 
as well as the interaction energies. In the following, we 
focus on these pairwise interactions, for two generic Burg- 
ers vectors. For configurations where the sum of Burgers 
vectors vanish (such as cylinders with periodic bound- 
ary conditions along the cylinder axis) the divergences 
discussed above will cancel. 

As noted in Ref. [9 , the Laplacian on a cylinder 
is equivalent to that in an infinite two-dimensional flat 
space together with the periodicity requirement. We ex- 
ploit this idea by first considering the stresses without 
the periodicity requirement. Up to a sign depending on 
the direction of the Burgers vector b — ±bx, the results 
for the infinite flat space are given by pQ: 



is different from a three-dimensional isotropic solid [18) . 
For a dislocation with b — ±by one has to take x — > y and 
y — > —x on the right hand side of the above equations. 

On a cylinder, a dislocation at (xo,yo) — (0,0) must 
be duplicated at intervals of W in the y direction, to 
respect the boundary conditions. Therefore the stresses 
at a point (x, y) generated by a a dislocation with b = bx 
at the origin are given by: 



Y> -Ab[y + kW][3x 2 + (y + kW) 2 ] 

1^^ [ x 2 + ( y + fcW)212 ' ( '' ) 



oo 

E 



Ab(y + kW)[x 2 - (y + kWf 
[x 2 + (y + kW) 2 ] 2 



(G) 



and 



x _ Abxjx 2 - (y + kW) 2 

Txy ~ fe ±"oo + (y + kw?f 



(7) 



The superscript on the stresses reminds us that these are 
the stresses created by a dislocation with b = bx. Similar 
sums represent the stress produced by an edge dislocation 
with b = iby. 

The force on another dislocation at (x,y) due to this 
stress will then be given by the Peach-Koehler force [4j 
HUES]: 



Fi — bko-jk^i 



(8) 



where €ij Z is the Levi-Civita tensor. Explicitly, we have 
for b = bx: 



Jlat _ 



= -Aby(3x 2 +y z )/r 4 



(2) 
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and 



= Aby(x 2 - y 2 )/r 4 



T flat _ 
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Abx(x 2 -y 2 )/r 4 



(3) 



(4) 



while for b = by: 



Fy b(T X X ■ 



(10) 



withr 2 = x 2 +y 2 and A = where Y = 4^(/!+A)/(2/z+ 
A) is the two-dimensional Young's modulus. Note that 
the functional form of the spatial dependence, being de- 
termined by geometry, is the same as for plane stresses 
around dislocations in three-dimensional solids. How- 
ever, the relationship between A and the elastic constants 



F X b(Tyy , 



Fy b(7 x y • 



(11) 



(12) 
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Summing the series 

The sums of Eqs. ^ through ^ can be evaluated 
using the Sommerfeld- Watson transformation [36]. To 
demonstrate this for the first sum, consider the function 
g(z) = cot(7rz), which has only simple poles of residue 
unity which lie on the x axis at integer values. The sum 
of Eq. ([2]) can then be written as the complex contour 
integral: 



substitution x 
to: 



y, y 



-x, which leads immediately 
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(19) 
(20) 



Application of the Somerfeld- Watson transformation to 
the remaining component of the stress tensor leads to: 



g(z)f(z)dz, 



(13) 



with: 



(y/W + z)[3(x/W) 2 + (y/W + z) 2 } 
J[ ' W[(x/W) 2 + (y/W + z) 2 } 2 



(14) 



Since f{z) ~ 1/z at large distances from the origin, we 
can deform the contour so that it captures only the poles 
of f(z) (note that the integral on the circle at infinity 
vanishes even though the decay is only ~ 1/z, due to the 
cot(7rz) term). Upon rewriting the function f (z) as: 



f(z) = -(Ab/W) 



(y/W + z)[3{x/W) 2 + (y/W + zf 



(z + (y + ix)/W) 2 (z + (y- ix)/W) 2 ' 

(15) 

we see that it has two poles of order 2. Summing the 
residues gives: 



•r; , = j> g(z)f(z)dz 



2W 2 
■nAb 

2W 



- (csc 2 (7r(y-ix)/W) -csc 2 (7r(y + ix)/W)) 
(cot(7r(y + ix)/W) + cot(7r(y - ix)/W)) . (16) 



In a similar fashion one obtains the other components of 
the stress produced by a dislocation with b = bx: 



iAbn 2 x 



(esc 2 (tt (y - ix) / W) - esc 2 (tt (y + ix) / W) ) 



2W 2 

— (cot(7r(y + ix)/W) + cot(7r(y - ix)/W)) , (17) 



Abir 2 x 
' 2W 2 



(csc 2 (7r(y - ix)/W) + csc 2 (?r(y + ix)/W)) . 

(18) 



To find the components of the stress tensor due to a 
dislocation with b = by, we use the previously mentioned 
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_ ^7 7r [ _2cotn [ 7r ( :z: ~ iy)/W] - 2coth[7r(x + iy)/W] 

7T— (csch[7r(x - iy)/W} 2 + csch[7r(a: + iy)/W] 2 )]. 

(21) 



Asymptotic forms 

It is natural to consider various limits for the stresses. 
For distances small compared to the cylinder's radius, it 
can be checked that the previous expressions all reduce to 
the fiat space results, as must be the case. However, for 
large separations along the x direction (i.e., the cylinder 
axis), the behavior is different. For a dislocation with 
b= bx, we find that for \x\ > W = 2ttR: 



a x xx a -47r 2 A6e- 2x|a;|/M/ (|a ; |/V7 2 )sin(27ry/W), (22) 



a x yy w 4^ 2 A6e- 2T l 2; l/ w/ (|a;|/iy 2 )sin(2^y/W), (23) 



a x xy » 47r 2 A6e- 2w l K l/ w (|a;|/^ 2 )cos(27ry/W), (24) 
while for a dislocation with b = by, we have: 



a y xx a 4Tr 2 Abe- 2 ^ xl/w (\x\/W 2 )cos(2Try/W), (25) 



a'l y « 2nA(b/W) ■ Sgn{x) 

- 47r 2 A6e- 2,r|a;|/M/ (|a;|/VF 2 )cos(2^y/W), (26) 



0% w ATT 2 Abe- 2 ^ xl/w (\x\/W 2 )sm(2iry/W), (27) 

where Sgn(x) = x/\x\. Thus, all components of the stress 
tensor decay exponentially, except for the circumferen- 
tial stress induced by a dislocation with b = by, which 
approaches to a constant value exponentially fast. This 
constant reflects the half-line of extra material that is in- 
serted throughout the long axis of the cylinder, leading 
to a long-ranged stress field, as evident in the first term 



ofEq. (26) 
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Energy considerations 



In this section we convert our results for the stresses to 
the elastic interaction energies for two dislocations on a 
cylinder, obtained by integrating the Peach-Koehler force 
of Eq. As discussed above, a generic Burgers vector 
can be decomposed into x and y components, and thus 
we consider three distinct scenarios: 

(a) Both dislocations have Burgers vectors in the ±x 
directions. 

(b) One dislocation with b = ±bx with another with 
b = ±by. 

(c) Both dislocations have Burgers vectors in the ±y 
directions. 

In flat space, the interaction energy of two edge dislo- 
cations with vectors b\ and 62, and with a relative sepa- 
ration of r = (x, y), r> b, is given by [5]: 



E(x,y) = -A[(b 1 -b 2 )log[r} 



(fti-r)(fr 2 -r) 



2E r 



(28) 

where b is the lattice spacing, and the effect of higher 
order terms in the gradient expansion of Eq. ([I]) is given 
by the core energy term 2E C [37]. Note that unless the 
Burgers vectors are equal and opposite, the total energy 
of the system will also include terms that diverge with 
the system size, as previously discussed. 

The derivatives of the interaction energy with respect 
to the coordinates yield the forces: for example, differ- 
entiating Eq. (28) with respect to x or y and using Eqs. 



@ and (TO) gives Eqs. Q and Up to a constant, we 
can obtain this interaction energy by integration of the 
Peach-Koehler force. 

Upon generalizing to the case of the cylinder, where 
we have already found an explicit formula for the forces, 
we can use it to derive the expression for the interac- 
tion energy for case (a) above. Integrating the force in 
the x direction, F z 
E(x,y)=Y(y) + C(x,y) with: 



ba% y , with respect to x leads to 



C(x,y) 
Ab 2 



Ab 2 



log[sinh(7r(x — iy)/W)\ + 



+ 1 ^-i7r(x/IU)csc(7ry/W)sinh(7rx/W)csch(7r(x — iy) /W) 

+ C.C., (29) 

with W = 2ttR. The derivative of C(x, y) with respect 
to y can be shown to be equal to F y — ba xx , implying 
that Y(y) = const. The constant can be found be de- 
manding that the expression reduces to that of flat space 
for W 3> x,y, see Eq. (28). Our final result for antipar- 



Ab 2 W 

E^i{x,y) = — log[— sinh(7r(a; - iy)/W)\ + 
Ab 2 

+ -^-i7r(a:/IU)csc(7ry/W)sinh(7rx/W)csch(7r(x — iy)/W) 

+ C.C. (30) 

We have suppressed, for simplicity, the large distance 
core energy contribution displayed in Eq. (28). The no- 
tation E±,~x denotes that this is the interaction energy 
of two antiparallel dislocations with Burgers vectors in 
the ±x directions. For x ^> R, we find that: 



Ex 



-ifarH/w C0S (27ry/W). 

(31) 

Fig. [5] shows the equal energy contours of this inter- 
action energy. Close to the origin, a cut parallel to the 
x-axis would give a graph with two minima, correspond- 
ing to the two dislocations with separation vector at a 45 
degrees angle to the i-axis - the double minima are ex- 
pected, since close to the origin we are not sensitive to the 
finite radius of the cylinder, and this is indeed the stable 
configuration of two dislocations in flat space, when climb 
processes (motion perpendicular to the Burgers vector) 
are prohibited [38] . The double minima structure at a 
fixed offset y = W/10 is shown in Fig. |6] However, at 
larger vertical separations the minima become shallower 
and shallower, until at a separation of W/4 along the cir- 
cumference, there is only a single maximum at x = 0. 
For fixed y, the minima obey d-Eg.-^.y) x a x _ q 



Upon equating Eq. ( 18 ) to zero, we find (aside from the 
solution x = 0): 



tan(?ry/W) = ±tanh(?rx/W), 



(32) 



which indeed has a solution at nonzero x provided y < 
W/A. 



Taking the y — > limit of Eq. ( 30 ) (which is ill-defined 
for y = 0) leads to: 



Ez-£(x,y-> 0) = 



r 2R 



Ab 2 log[— sinh(7rx/W)] - — coth(7rx/W) 



W 



(33) 



This formula is equivalent to Eq. (2.1b) in Ref. |38j . 
describing the interaction energy of two antiparallel grain 
boundaries in the infinite plane. The connection between 
these two systems will be elucidated in section [J 

In a similar fashion we can find the interaction energy 
for case (b) above, of an edge dislocation with b = bx with 
another with b = by. The forces on the latter are given 



allel Burgers vectors along the cylinder axis is thus: 



by F x 
find: 
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-o%y By integrating the stresses we 





FIG. 5. Equipotential contours for a dislocation on a cylin- 
der at position (x, y) with b = — bx, interacting with another 
dislocation at the origin with b = bx. Energy is measured in 
units of Ab 2 . Note the much lower energies than in Fig. [8] 
where the Burgers vector are rotated by 90°. 
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FIG. 7. Equipotential contours for a dislocation at (x, y) on 
a cylinder with b = by, interacting with another dislocation 
at the origin with b = bx. Energy is measured in units of 
Ab 2 . Note that the large distance core energy contribution, 
analogous to that appearing in Eq. (28 1, is undefined with- 



out additional dislocations on the cylinder, since the sum of 
the Burgers vectors should vanish. Nevertheless, this expres- 
sion is useful for charge-neutral dislocation configurations on 
a cylinder, whose interactions can be decomposed into pair- 
wise interactions involving Ex,-i, j 
(|3TL(|36l) and JMb. 



0. For x, y < W, Eq. 



For y 



0, we find that E xy 



( 34 ) reduces to the fiat expression: 
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Ab 2 xy 
x 2 + y 2 



Ab 2 
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FIG. 6. Slice along the energy contours of Fig. [5] showing 
interaction energy as a function of x, for several ratios of y/W . 
Note the double minima, which become a single maximum for 
V > W/4 



Ab 2 n^sm[2ny/W} 
coa[2wy/W] ~ cosh[27ra/VF] ' 1 ' 



The equal energy contours of this interaction energy are 
shown in Fig. [7] For x ^> R the expression reduces to: 



w 2ttA6 2 — sm[2%y/W 



-2tt\x\/W 



w 



Since \c + ~\ > 2 for any value of c, there is clearly 
no divergence of the energy as the two dislocations come 
together, with a minimum energy of — 4|- obtained for 
x = —y. The fmiteness of the energy as the dislocations 
as (x, y) — > (0, 0) reflects the vanishing of the logarithmic 
divergence in Eq. ( 28 ) when bi and & 2 are perpendicular. 



As mentioned before, a dislocation configuration where 
the sum of Burgers vectors does not vanish leads to ad- 
ditional terms in the expression of the total energy di- 
verging with the system size. The above equations only 
reflect the pairwise interaction terms. 

Finally, in case (c), the force on an edge dislocation 
with b = —by induced by another dislocation with b = by 



is given by F x 



yy 



F„ 



<j v 



Upon integration we 




Stresses and energetics for Burgers vectors at 
arbitrary inclination angle 

A Burgers vector b that makes an angle 9 with the 
x axis, can always be decomposed into Cartesian com- 
ponents parallel and perpendicular to the cylinder axis. 
Since we assume linear elasticity, the relevant stress fields 
follow from the superposition of the solutions obtained 
previously: 



(39) 



appear in Eqs. ( 16 21 1 



where afj and cr^ 

We now determine the interaction energy of this dis- 
location with another dislocation, whose Burgers vector 
forms an angle a with the x axis (on a triangular lat- 
tice, if the Burgers vectors have their minimum allowed 
lengths, the difference of the two angles 9 and a will be 
a multiple of 7r/3). The force on this dislocation is then: 



FIG. 8. Equipotential contours for a dislocation at (x, y) on 
a cylinder with b — —by, interacting with another dislocation 
at the origin with b = by. Energy is measured in units of Ab 2 . 
Note that for x 2> R — W/2n the energy contours become 
parallel to the y axis; so that the forces are primarily in the 
x direction. This trend, and the associated linear potential, 
can be seen from Eqs. |TT|, (p2|,j26| and p7||. 



find 



Fr 



b[a xy cos 9 cos a 



a yy cos o sin a 



a v yy sin 9 sin a] . (40) 



F, 



y — b[—o% x cos 9 cos a 



a xy cos V sin a 



Upon integration we find that: 



a v xy sin 9 sin a] . (41) 



Ab 2 W 

E v-y( x ,y) = ^rlog[^-sinh[7r(a: - iy)/W\\ 
Ab 2 

— in (x/W) esc [ny / W] sinh [nx/ W] csch[-7r (x — iy)/W] + 

C.C., (36) 

where e is the base of the natural logarithm. The equal 
energy contours are shown in Fig. [8j The analytic form of 
Eq. (36 1 is similar (but not identical!) to that for E x _ x , 



see Eq. p0|). As before, we can check that for x,y <C R 

(this correspondence requires 
For x R we 



this reduces to Eq. p8| ) 
the extra factor 2/e in the logarithm), 
now find, in contrast to Eq. (31 ), a linear potential, 



E 



y,-y 



2-kAK 2 
W 



x\. 



(37) 



En 



E x 



. cos f cos a - 



- E„ 



- Ea, 



, sin f sin a. 



(42) 

For the phyllotaxis problem discussed in section]] the 
energy landscape associated with dislocation pairs with 
antiparallcl Burgers vectors is of particular interest; these 
can nucleate locally and then unbind, thus modifying the 
geometry of the lattice. Fig. [9] shows an example of the 
energy equipotential contours for 9 = 7r/6 and a — 7ix /6. 



Structure of the energy landscape 

In this section we describe in more detail the structure 
of the interaction energy landscape for two dislocations 
interacting on the surface of a cylinder at zero tempera- 
ture. 



In the limit y — > 0, similarly to the case of E$ t -$(x, 0), 
we find that: 



Burgers vectors are in the ±x directions 



Ey,_y(x,0) = Ab' 



W 

log — -sinh(7rx/W)] 
e7rb 



W 



3th(7rx/W) 
(38) 



In the case of dislocations with antiparallel Burgers 
vectors in the ±i directions, the interaction energy is 



given by Eq. (30). From the energy contours (see Fig. 
[5]), we can see that the configuration where the two dis- 
locations are located at the same x coordinate but on 
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x/2ttR 



FIG. 9. Equipotential contours for a dislocation on a cylinder 
with a Burgers vector b — &(-^ai + forming an angle 
of 7r/6 with respect to the i-axis, interacting with another 
dislocation with an opposite Burgers vector. The interaction 
energy is given by Eq. (42 1. Note the two skewed saddle 



points. Energy is measured in units of Ab 



opposite sides of the cylinder, the anti-podal point, is 
a maximum of the energy. If the two dislocations have 
the same sign, there is a stable minimum riding on top 
of an infinite energy due to the lack of overall charge 
neutrality. To see this, one can expand the force on a 
dislocation offset slightly from the origin by an amount 
(x, y), exerted by another dislocation at the anti-podal 
point. This gives, for antiparallel Burgers vectors, to 
lowest non-vanishing order: 



F x (x,y) = Kx, 
F y (x,y) = Ky, 



(43) 



with K = 



> 0, so the minimum is unstable. For 
parallel Burgers vectors, the signs are reversed. 

We conclude that two dislocations with antiparallel 
Burgers vectors on opposite sides of the cylinder can min- 
imize their energy by annihilating, while two dislocations 
of the same sign will remain at the anti-podal point to 
minimize the repulsive elastic interaction energy. 



Burgers vectors are in the x and y directions 

In this case the anti-podal point is found to be a saddle 
point, as can be seen near the center of Fig. [7] As before, 
we suppress an infinite energy due to the lack of charge 
neutrality. 



Burgers vectors are in the ±y directions 

Upon repeating the same analysis for two dislocations 
with Burgers vectors in the ±y directions, we find that, 
again, the only extremum point is at the anti-podal point 
(0, W/2), which Fig. [8] shows to be a saddle point. 



NUMERICAL TEST ON A TRIANGULAR 
LATTICE 

As discussed in the Introduction, in soft matter physics 
a triangular lattice of colloidal particles wrapped around 
a cylinder may be realizable. If one of the principal axes 
of the lattice is oriented along the circumferential direc- 
tion, the minimal Burgers vector of a dislocation can take 

on the following values: b — ±by, b = ±6 (^^x ± \\ 

To test our continuum limit predictions numerically, 
consider the stresses associated with an isolated disloca- 



tion at the origin with Burgers vector b = b( 



VI; 



hv) 



(i.e. — — 7t/6) in a triangular lattice of masses and har- 
monic springs which are relaxed to their minimum energy 
configuration by a conjugate gradient method. To insure 
a charge-neutral configuration, an additional dislocation 
with antiparallel Burgers vector was created at the nega- 
tive x end of the cylinder. The strains and stresses at any 
point can be calculated from the shift in position of the 
surrounding points relative to the perfect lattice. In Fig. 
(10), wc display a xx and a yy for positive a; at a constant 
y = W/4 for cylinders of various sizes. We obtain good 
agreement with our continuum results even for relatively 
small system sizes, which can be realized in colloidal ex- 
periments. This calculation illustrates the applicability 
of our results for Burgers vectors rotated away from the 
x-axis, which lead to long range strain fields along the 
cylinder axis. 

We note from Fig. ( f ) that the asymptotic value of 
a yy at large positive x in the simulations is zero rather 
than ayy' e (x — > oo) = — nA/W predicted from Eq. (39 1 . 



This is due to the contribution of the additional disloca- 
tion positioned at the negative x end of the cylinder: 
from Eqs. (23), (22), (24) we see that while it has an ex- 



ponentially small effect on the other components of the 
stress tensor, it creates a constant circumferential stress 
a yy , thus shifting the values for the circumferential stress 
by a constant. At x 3> W, the stress created by the dis- 
location at the origin is also approximately constant, and 
since the two dislocations are antiparallel their contribu- 
tions are equal but of opposite sign - which explains why 

a yy = a yy^ + a vv^ — > at large positive x. 

There is some evidence for helical motion of tracer par- 
ticles on the cell walls of elongating bacteria such as Es- 
cherichia coli |39j . suggesting that the lattice shown in 
Fig. |3]may be slightly skewed. Hence, it is of some inter- 
est to find the asymptotic form as x — > oo of the stress for 
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FIG. 10. Stress components a xx (left) and a yy (right) as a function of a; at a constant offset around the cylinder y = W/4, for a 
dislocation with 6 = b(^x — |y) at the origin. To account for a charge- neutral dislocation configuration, another dislocation 



with an antiparallel Burgers vector is placed along the negative x end of the cylinder, far enough from the origin to create 
a constant o yy stress [see Eq. |26l] and a negligible contribution to the other components of the stress tensor. The symbols 
are for simulations of cylinders with W increasing from 86 to 406 where 6 is the equilibrium lattice constant of the simulation 
lattice. In all cases, the length of the cylinder is 8.66VF and the dislocation is situated at the center of the cylinder along the 
long axis to minimize edge effects. The solid line is the theoretical prediction from Eq. (39 1 with 6 — — 7r/6, where for the right 
figure one has to add the constant a yy stress mentioned above. Except for one or two points closest to the dislocation, where 
discrete effects become important, the simulations agree with the results of continuum elasticity even for relatively small system 
sizes (W > 166). For negative values of x, the convergence to the continuum limit results is slower, and there are corrections to 
the stresses which scale as ~ b/W. These will be discussed in more detail in future work. 



a dislocation with a Burgers vector in a direction forming 



a small angle 9 with the x axis. We use Eqs. (26) and 



(391 to find: 



<J yy ~e- 2irAb/W ■ Sg(x) 

+4ir 2 Abe- 2 ^ xl/w \x\/W 2 [sm(2iry/W) - 0cos(27ry/W)]. 

(44) 

Thus, for slightly tilted Burgers vectors there is a 
length scale I* along x at which the constant contribution 
to the stress dominates over the exponentially decaying 
one, which is readily found by equating the two terms: 



[2ire- 27rr/w r/W}sin(2iry/W) = 0. 



(45) 



Provided y is not too small, with W = 2wR, we find 
that to leading order: 



I* ~ Rlog(0). 



(46) 



ANALOGY WITH GRAIN BOUNDARIES 

We now describe a correspondence between disloca- 
tions on a cylinder and grain boundaries, which allows 
checks of some of our calculations 00] . 



Suppose that the lattice orientation on the cylinder al- 
lows two dislocations of opposite Burgers vectors, ±bx, 
initially very close together. Let us now glide one of the 
dislocations away from the other, until they reach a cer- 
tain separation A along the x, which coincides with the 
axis of cylindrical symmetry. Every time the dislocation 
glides by one lattice spacing, we connect an initially cir- 
cumferential row with an adjacent one to make a spiral. 
Thus, at any fixed separation, the finite part of the cylin- 
der between the two dislocations is converted to a helical 
structure, which can be thought of as a local rotation of 
the crystalline lattice (see Fig. |4j). Since a shift of b is 
associated with each rotation, the pitch of the helix is 
6 = b/W . Therefore, we can view the circumferential 
row passing through each of the dislocations as a grain 
boundary, and the whole structure as a polycrystalline 
material with three regions, the middle one tilted by 9. 
Glide separation and the equivalent grain boundary con- 
figuration are also illustrated in Fig. [II] as well as in 
the Supplementary Video SI, showing numerical results 
for a system of masses and springs, where the minimal 
energy configuration of the system at each instance in 
time is obtained by a conjugate gradient minimization 
procedure. 

By this analogy, the interaction energy E(x, y) of two 
dislocations on the cylinder separated by a displacement 
(x,y), corresponds to the energy (per section of length 
W) of two infinite, straight, grain boundaries a distance 
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FIG. 11. Dislocations on a cylinder are analogous to grain boundaries. The images show configurations obtained numerically 
for a hexagonal lattice of masses and harmonic springs whose energy is minimized as the dislocation pair is separated by 
a glide. On this triangular lattice, a dislocation corresponds to a point with five neighbors adjacent to a point with seven 
neighbors. Starting with two nearby dislocations, (a), we glide one of the dislocations to the right in steps of the lattice spacing 
[Supplementary Video SI; figures (b-d) show intermediate snapshots during the glide]. The finite section of the cylinder between 
the two dislocations displays a helical structure, similar to the model with magnetic beads shown in Fig. [4] The analogy with 
a flat polycrystalline material is illustrated by unrolling the cylindrical crystal in (d) and duplicating it several times in the 
y (circumferential) direction to replicate the periodic boundary conditions (e). The resulting infinite columns of dislocations 
(dotted lines) define two grain boundaries separating the middle region, rotated by an angle 6 = b/W, from sections with the 
original crystal orientation on either side. 



x apart and offset by an amount y along their lengths, 
where the distance between adjacent dislocations in each 
grain boundary is W. As is well known in the theory of 
small angle grain boundaries [H [38] , the associated grain 
boundary angle is = b/W. For example, the interaction 
energy of two parallel grain boundaries with Burgers vec- 
tors in the ±i directions and separated by a distance x, 
mathematically equivalent to the situation in Fig. 11 'c), 
is given by Eq. (33) ; This is in agreement with Eq. 
(2.1b) in Ref. [38]. 



ENERGY LANDSCAPE IN THE PRESENCE OF 
A CONSTANT AZIMUTHAL FORCE 

Consider two dislocations with antiparallel Burgers 
vectors ±bx on a cylinder, a small distance apart. For 
simplicity, we pin one dislocation at the origin. Let us 
assume that there is an additional constant force acting 
on it, due, say, to the Peach-Koehler force created by an 
external stress a xx [see Eq. (10)]. 



The energy landscape now resembles a tilted wash- 
board potential with period W and is given by: 



E x-x( X 'V) 

Ab 2 



Ab 2 W 

— - log — ; sinh(7r(a: - iy) W)\ + 

2 7TD 



wr(x/W)csc(7ry/W)sinh(7rx/W)csch(7r(x — iy)/W) 

+ C.C. - Gy, (47) 

where the constant azimuthal force G oc a xx . For climb 
dynamics, there can also be a contribution due to the 
chemical potential associated with adding new material 
[S] . The dynamics of the dislocation pair in this potential 
will depend on the glide and climb mobilities to be intro- 
duced below. First, however, we discuss how the energy 
landscape changes due to the field G. Clearly, the mini- 
mum at the origin remains if the field G is not too large. 
Now, however, there are important new saddle points in 
the energy landscape. Let us first locate these saddles in 
the case of flat space. Here, the potential is given by Eq. 
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(28), with an additional term —Gy. A necessary condi- 



tion for the existence of a saddle is the vanishing of ^ 
and Upon differentiation of _~{x,y) we obtain: 



^saddle Vsaddle — 



G 



(48) 



This extremum must be a saddle point since (a) for a 
given y coordinate there is a minimum energy at x — y 



[see discussion preceding Eq. (32)], and (b) this cannot 
be a minimum in the y direction since y — > and y — > oo 
give us a negative, diverging energy. 

Let us now consider the problem on a cylinder. The 
only scales that may enter the problem are Ab 2 , W and 
G; b cannot enter explicitly since we are considering the 
continuum limit. Only one dimensionless parameter can 
be formed out of these, namely: D = GW/Ab 2 . There- 
fore the y position of the saddle (which as we show below 
exists for any value of G), can be written as: 




FIG. 12. 
with 



aw 



Equipotential contours for a dislocation pair 
1, with W — 2nR. Note that the force G caused 
by a macroscopic stress a xx has pulled the maximum away 
from y = W/2. Note also the two saddle points (denoted by 
the letter S) near y « W/4. 



Ab 2 

Vsaddle = ~^-f(D), 



(49) 



where f(D) is a function of the dimensionless parameter 
D. The case of flat space, Eq. (48), corresponds to D — » 



oo, showing that f(D) must asymptotically approach the 
value one. We now proceed to discuss the general form of 
f(D). The condition ^ does not depend on the field in 
the y direction, and thus, Eq. (32) still holds. Therefore, 



we search for the saddles whose coordinates have the form 
(z*{v)>y)> where: 

W 

x*(y) = ±— arctanh(tan(7ry/W)), (50) 

with < y < W/4. 

Within this parametrization, we still have to satisfy 
= 0. This condition implies that F y — G — ba% x = 0, 
with given by Eq. ( 16 ). For small values of y, a van- 
ishing derivative in the x direction would yield x*(y) = y, 



as can also be seen from Eq. (50). For such a point, we 



know that the two dislocations have a large force in the 
y direction pulling them together, which can be seen di- 
rectly from Eq. ^ (since in this regime we are not sen- 
sitive to the cylindrical geometry). On the other hand, 



for y very close to W/4, Eq. (50) yields an x coordinate 



which diverges. Hence, the stress in Eq. ( 16 ) must van- 
ish, since the force between two dislocations falls off as 
a power-law. This argument implies that the y compo- 
nent of the force at such a point will be approximately 
G, and in particular, it will be positive. Therefore, there 
is an intermediate value of y for which both F x and F y 
vanish. We conclude that a saddle exists for any value of 
G! (In fact, two saddles, since we have a reflection sym- 
metry around the y axis.) Equipotential contours where 



for values = 1 and a much stronger force given by 
= 10, respectively. For large values of G, the saddles 
coincide with those obtained for flat space. However, for 
small G the saddles will be close to the line y = W/4 and 
their x values will diverge. Fig. ( 14 ) shows this scaling, 
expressed in Eq. (49), and the asymptotic limits: for 
GW/Ab 2 — > oo, the rescaled y coordinate of the saddle 
approaches unity, corresponding to the flat space re- 

Hp-, corresponding 



Ab 2 

suit. For the opposite case, jjp- 



to the straight line f(D) = D/4. To determine the form 
of the divergence of the x position of the saddle for small 



G, note from Eq. ( 22 ) that at large values of x the force 



in the y direction decays exponentially for any value of 
V, as: 

F y a [4^ 2 A6 2 e- 2,ra/w Viy 2 ]sin(27ry/W). (5 1 ) 

Upon equating this force to G, we find for small G a 
saddle at x* ~ -^( ga^p; ), a location 
that diverges approximately logarithmically as G — > 0. 

An important feature of the saddle when considering 
the thermal excitation of dislocation pairs is its energy 
relative to, say, the interaction energy of the two dislo- 
cations when they are one lattice spacing apart. This 



energy difference is plotted in Fig. [15j For flat space, it 
is easy to see that the energy barrier diverges logarith- 
mically as G — > 0, and is given by: 



U(G) ps Ab 2 log(A6/G), 



(52) 



the saddles can be seen are shown in Figs. ( 12 ) and ( 13 1, 



where the energy is measured relative to that of a dislo- 
cation pair separated by a lattice constant b. 
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FIG. 13. Equipotential contours for a dislocation pair with 
j^gz — 10. Note the maximum and two saddle points near 
the line y = W/W. 



Eq. (52 1 is approximately correct on the cylinder as 
well, provided the y coordinate of the saddle is much 
less than W. For small values of G, more careful anal- 
ysis is required. Let us parametrize the saddle point by 
(x*(y),y). From Eq. (18) we infer that at large values 



of x the force in the x direction decays exponentially for 
any value of y, as: 

F x « [An 2 Ab 2 e- 2 * x/w x/W 2 ]cos{2Try/W). (53) 

Hence, the potential at infinity does not diverge, leading 
to a finite energy barrier for dislocation unbinding at any 
value of G, as shown in Fig. (15). We can adapt the 



scaling analysis for the y coordinate of the saddle also 
to the energy. Consider the interaction energy E(x, y) in 
the absence of G. Although the lattice spacing b enters 
the expression for the interaction energy explicity, it can 
be easily eliminated: from Eq. (31 ) we see that 




Cylinder 

Flat, y = G/Ab 2 

y=w/4 



aw 



FIG. 14. The rescaled y coordinate of the saddle point ^=jj 
is shown as a function of the dimensionless parameter §K . 
The horizontal line is the result in flat space, where y = 



- Cylinder 

-Flat, U = Ab 2 log(Ab/G) 




FIG. 15. The saddle point excitation energy U, given by Eq. 
(31 1, is rescaled and plotted as a function of ^*r. In the flat 
space limit the scaling function is given by r/(D) = — log(D), 
corresponding to U = Ab 2 log(Ab/G). For the opposite case, 
when G is small and the cylindrical geometry is important, 
the energy U saturates at a constant value even as G — ► 0, 
const + log(f ). 



E(x, y) = Ab 2 log(W/6) + Ab 2 h[x/W, y/W], (54) 

which can also be deduced from dimensional consid- 
erations combined with the flat space limit. The energy 
U(G) is thus given by: 



U(G) = E(x saddle, y saddle) - E{\f\ = b)~Gx s 



uldlc 



(55) 



Previously we found that the y coordinate of the saddle 
point scales as y sad die = ^ff{D), see Eq. (49 1. The 
same argument can be repeated for the x coordinate of 
the saddle, giving x sa ddie — ^%-gi.D), where the function 
g also depends on the dimensionless parameter D. Plug- 
ging these results into Eq. ( 55 1 and using Eq. ( 54 ) we 
find that: 



U{G) = Ab 2 r,(D) - Ab 2 log(b/W), 



with ?y(Z?) = h[2^p-, ^ffi ]—g(D) a function of the dimen- 
sionless parameter D. Rearranging the above equation, 
we find a dimensionless scaling form: 



(56) 



U(G) - Ab 2 log(VF/6) 
IP 



V(D), 



(57) 



For flat space, we have U(g) = Ab 2 \og(Ab/G), showing 
that for large values rj(D) w — log(Z?). The above form 
of scaling is illustrated in Fig. [T5j 
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THERMAL NUCLEATION OF DISLOCATION 
PAIRS: A TWO-DIMENSIONAL ESCAPE OVER 
A BARRIER PROBLEM 



With colloidal particle arrays or bacterial cell walls at 
relatively high temperatures in mind, we shall now cal- 
culate the thermal nucleation rate of dislocations: in the 
presence of a uniform external field denoted G, two dis- 
locations of opposite Burgers vectors can overcome the 
energy barrier due to their attractive interaction and un- 
bind due to thermal fluctuations. We shall first treat 
this process in flat space, and then compare to the re- 
sults in cylindrical geometry, which differ significantly 
when G is small. Our analysis will be done for the case 
of dislocations with Burgers vectors ±6x, where the ex- 
ternal force driving the nucleation is in the y direction, 
perpendicular to the Burgers vectors. For simplicity, we 
assume isotropic elastic constants. The same strategy 
can be adapted to other scenarios. We then discuss the 
nucleation of dislocations along the cylinder axis due to 
a twisting stress a xy . 



Nucleation of dislocations in flat space 

We consider first the thermal escape-over-a-barrier 
problem for two dislocations in flat space (similar results 
should hold for a cylinder with a large enough field G 
such that the saddle distance from the origin is much 
smaller than W). It is convenient to assume that one of 
the dislocations is mobile, while the other is pinned at the 
origin. Alternatively, we can consider the motion relative 
to the fixed center-of-mass of the dislocation pair. For a 
similar analysis of thermal activation of vortex pairs in 
superfluid helium films, see Ref. [51]. Throughout this 
section, we choose units such that the Boltzmann con- 
stant ks = 1- 

We shall assume over-damped dynamics for the mobile 
dislocation, described by an anisotropic diffusion tensor 
D [42]. We analyze the scenario where the dislocations 
can glide and climb. This choice is motivated by the bac- 
terial elongation problem discussed in section [J where it 
is the climb of dislocations that drives the growth pro- 
cess. The motion of dislocations for colloidal arrays on a 
cylinder could also have a climb component, provided the 
particles can jump on and off the cylinder, thus providing 
an external source of vacancies and interstitials. Notice 
that allowing climb mobility is very different from the sit- 
uation in Ref. |38j , where dislocations can only glide and 
one has a one-dimensional escape-over-a-barrier problem, 
rather than the two-dimensional problem that we study 
here. 

The probability n(r, t) of finding the dislocations with 
separation r is given by the continuity equation: 



dn 

— + dwj = 0, 
dt 

with a probability current j(f,t): 

j = finF(r) — DVn. 



(58) 



(59) 



The first term is given by F{r) — — VC/(r), where the 
potential energy U(r) is the sum of the interaction of the 
two dislocations as given by Eq. (31) and the potential 
— Gy associated with the external driving force. Suffi- 
ciently close to equilibrium, the Einstein relation relates 
the diffusion tensor to the mobility tensor: D = T^t, 
allowing us to recast the Fokker-Planck equation as (set- 
ting k B = 1): 



dn 



-V[-De- c/ / T V( 



lie 



(60) 



We now exploit Langer's formalism [T^] to determine the 
structure of the probability currents that describe the 
nucleation of dislocation pairs. 

For small enough escape rates (i.e., a low enough tem- 
perature) , the process will be dominated by currents flow- 
ing near the saddles points of the potential energy land- 
scape. In the case of a dislocation pair in flat space 
there are two symmetric saddles; the total escape rate 
will be double the probability current flowing through 
one of them. 

The escape rate through the saddle includes the Hes- 
sian (matrix of second derivatives) around it. We denote 
the Hessian in the vicinity of the saddle by H sa ddie, and 
the Hessian in the vicinity of the minimum by H m i n (for 
dislocation pairs, the latter requires a short distance cut- 
off of order the lattice spacing b). The probability density 
around the minimum follows a Boltzmann distribution, 
since thermalization in that region is a fast process com- 
pared to the time scales associated with crossing the high 
barriers. It is useful to define the transition matrix A as: 



DH 



saddle i 



(61) 



and denote its two eigenvalues by Ai > and A 2 < 0. 
The escape rate T is then given by [J_3] : 



r 



Ai / det[H mill } 



27T V I det [H sadd ie 



mini e~^ s 



e/T 



(62) 



We shall assume that D is anisotropic but is diagonal in 
a basis where one of the eigenvectors coincides with the 
direction of the Burgers vector, say, x. In this basis: 



D 



Dg 

D c 



(63) 
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where D c and D g are the glide and climb diffusion co- 
efficients. As mentioned before, for bacterial cell wall 
growth dislocation motion is predominantly via climb, 
i.e., D c ^> D g . For the case of colloids on a cylinder, we 
would typically have glide dynamics, i.e., D g 3> D c . 

Upon finding the positions of the two symmetric saddle 
points and calculating their Hessians, we obtain: 



H 



saddle 



G 2 
2Ab 2 



1 -1 
-1 -1 



It follows that \/|<2et[H sa ddie]| = 



V2Ab 2 



(64) 



and Ai = 



G 

. \b 2 



d 9 y 



8D C D 



Eq. (62 1 then leads to an escape rate given by: 



r oc 



Dg -D c 



with U Rl 



ddl e 



(Dc - D g y 



^ 2 l0g(^2) + C, 



8D P D B 



-U saddu /T 



(65) 

according to Eq. 



(31), with C a constant depending on the non-universal 



details associated with dei[H m in]> 

The proportionality constant for the rate in Eq. ( 65 ) 



will depend on the details of the dislocation nucleation 
at the origin, which should be independent of G. We 
conclude that the escape rate goes to zero as a power- 
law as G —> 0: 




FIG. 16. A mobile dislocation, initially at point A, interact- 
ing with a pinned dislocation with antiparallel Burgers vector 
at the origin O. The first dislocation is placed just beyond 
the saddle point, at (x,y) = (^j, j^), for — 10 (for a flat 
space, this would be the exact position of the saddle, however 
on a cylinder this point lies beyond it). The temperature is 
close to zero, although we assume motion over the Peierls po- 
tential [3] is possible. Here, W = 2-kR = 1, and the isotropic 
mobility tensor components are fi c = fi g = 1. After complet- 
ing one loop, the dislocation is "recaptured" and annihilates 
with the dislocation at the origin. 



r a 



G_ 



(66) 



. cos(27ry/W) sin(2^y/W) \ 
1 'l^Try/W) -cos(27ry/W) i " 



sm 



(69) 



Nucleation of dislocations on a cylinder 

We now solve the escape problem on a cylinder. The 
result for flat space should hold for large enough G so that 

A h 2 

^j=r- <C W. When this inequality is reversed, however, the 
saddles will occur (as discussed above) at a point with 
large x and y s» W/A, with U sa ddie approaching a finite, 
maximal value Ubarrier, independent of G (see Fig. 15). 
Con sider the form of H sa ddie in this regime. From Eqs. 



vTE 



(|24|) and \22\ we know the asymptotic form of ^ and 



dx 



Q y for large x. With this information we can calculate 
the Hessian, 



H 



saddle 




dF x 

6% 
dy 



(67) 



A straightforward calculation yields: 



H sad dle = [87r 3 Ab 2 27Tx/W 3 e'^^ w }M( y ), 



(68) 



where the matrix M is given by the y coordinate of the 
saddle: 



Since y « W/4, we have 

M RS 



1 

1 



(70) 



Using Eq. (62) we therefore obtain that: 



T cx y/D c D s 



sr(G)/T 



(71) 



Thus, as the field G becomes weaker, and the cylindri- 
cal geometry becomes more important, the escape rate 
decreases but saturates at a finite value! Note, however, 
that the escape problem on a cylinder does not fit pre- 
cisely to the conditions of the Kramers problem: one of 
the usual assumptions is that once the particle escapes, 
it can never be recaptured. Here, however, the opposite 
is true: the two dislocations can never truly escape from 
each other, due to the periodic boundary conditions. At 
low temperatures, the mobile dislocation is very likely to 
be "recaptured" and annihilated by its antiparallel part- 
ner after completing a single revolution, as depicted in 
Fig. |16| At finite temperatures, more revolutions are 
possible before recapture. 
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Nucleation of gliding dislocations along the cylinder 

axis 

Motivated by the bacterial elongation problem, in the 
previous example the dislocations had a finite climb mo- 
bility. For a crystal of colloids, however, the motion of 
dislocations is predominately via their glide, as previ- 
ously mentioned. With this motivation in mind, we set 
D c = in Eq. ( 63 ) and consider the unbinding of a pair 
of dislocation with b = ±x on a triangular lattice, where 
an external a xy stress tensor causes the dislocations to 
glide away from each other. Note that in spite of the pe- 
riodicity condition induced by the cylindrical geometry, 
this component of the stress tensor is not quantized, as 
explained in Appendix B. 

The o~ xy stress creates a force along the x axis, which 
is the glide direction. In this case, we have a one- 
dimensional escape-over-a-barrier problem, similar to 
that of Ref. [35], where the escape problem is solved 
in flat space. The potential in our case is different due 
to the cylindrical geometry, and is given by: 

Ab 2 211 

U(x) = -^-log[— sinh(7r(x — iy)/W)\ + 



Ab 2 



wr(x/W)csc(7ry/W)sinh(7rx/W)csch(7r(x — iy)/W) 



C.C 



ba xy x, (72) 



similar to the nucleation problem discussed above with 
G = ba xy . 

For large enough a xyi the maximum of the potential 
occurs for x — Abj a xy <C W, in which case the cylindri- 
cal geometry plays no role. Let us consider the opposite 
limit, namely, Wa xy <C Ab. In this case the maximum 
will occur for x max ^> W, where the potential decays 
exponentially. In infinite flat space, the barrier for nu- 
cleation would diverge logarithmically as o~ xy — > 0, as 
discussed in Ref. (38] . Here, on the other hand, the 
interactions are bounded, and thus the nucleation bar- 
rier approaches a constant U max . Upon using the one- 



dimensional version of Eq. (|62|) , we find that the prefac 



tor of the Arrhenius exponential term e u ' 
portional to \JU"{x max ), which vanishes as o~ xy 



is pro- 
-> and 

the maxima becomes shallower. Using Eq. (31 1 we find 
that that in this case the prefactor vanishes as ^/o~ xyi up 
to logarithmic corrections, so that: 



r oc 



(73) 



The above calculation represents a special case, since 
two antiparallel dislocations with Burgers vectors making 
a generic, finite angle with the x direction (say, and 
0— 7r), will have an interaction that provides a term linear 
in x for x » W. Using Eq. (37) and Eqs. (11) we find 



that the component of the interaction force in the glide 
direction is: 



F mt = 2irAb 2 cos(0) sm(0)/W 



(74) 



A a xy twist stress imposed on the boundaries of the cylin- 
der results in a force which has a component in the glide 
direction. Using Eqs. ([9| and (12) we find that: 



F g = F x cos(0) + F y sin(0) = ba xy [cos 2 (0) - sin 2 (0)]. 



(75) 



Thus, there will be a critical a* xy for which the inter- 
action term is compensated by the external stress, i.e., 

9 



2nAb 
W 



tan(20). 



(76) 



At very low temperatures, the dislocations would no 
longer be able to surmount the Peierls potential barriers. 
Then, the two dislocations will glide away from each only 
for a much larger stress, o~ xy ~ A, for which the force due 
to the external stress overcomes the periodic Peierls po- 
tential At finite temperature and a xy > a* y , we ob- 
tain a similar escape-of-a-barrier problem we had in the 
previous case, which can be solved in a similar manner. 
In this case, however, the glide direction forms an angle 
with the x-axis, and once the dislocations unbind they 
make a helical trajectory along to the glide direction. 
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The Airy stress function for a cylinder 

For some applications, it is convenient to describe the 
elastic stresses in terms of the Airy stress function x, de- 
fined by o~ij(r) = ei m ej n d m d n x(r) , which is a solution of 



the bi-harmonic equation [4J [18] . For a set of dislocations 
{b a } located at positions {r a } with zero external stress, 
the Airy function satisfies [43] : 



V 2 X (r) =4arA^2e ii b ct id j 5(f- 



(77) 



We can obtain x( x > V) from the strains, via 



dx 2 



'VV> dxdy 



-o~ X y Since for a dislocation with 
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b = bx, ba x xy = 
-b^-, we find: 



dxdy 



and 



dy 



Upon neglecting uniform translations, a general affine de- 
formation of this lattice can be written as: 



^- = -{l/b)E £ ,. x {x,y). 
Upon integration this leads to: 



Ab 

{irxy — a;iyarctan[cot(7ry/W)tanh(7rx/W)] 



(78) 



2W 

+(ix + y)Wlog 
Aby 



I _ e 27r(iy- X )/w' 



} 



(y + ix)log[sinh[7r(x - iy)/W] 

MW (~P 1 T 

zPolyLo 

4"7T V 



2 e -27r(x-iy)/W 



-c.c. 



(80) 



Although one could add an undetermined function 
G(x) to this expansion, the derivatives of x x with G(x) = 
satisfy the above equations, and thus Eq. (79) is the 
desired solution. 

Similarly, for a dislocation with b = by, we find: 



dx 



{l/b)Ey^y( X ,y), 



which upon integration leads to: 



x y {x,y) 

—iAby 



{Log 



^-2ir(x-iy)/W 



-Log(iSinhbr(:r + iy)/W])} 
Abe~™y/ W {~i + Tan[ny/W}) 



{2/Log 



4Sec[7ry/lU] 
{-ix + y)Log[Sinh[7r(a; + iy)/W]}} + C.C 



1 - e 



-2n(x-iy)/W 



(82) 



Although adding an undetermined function.// (y) is 
again possible, one can check that for H(y) = we have 



required. 



VV> dy 2 ( 2 '' y) ~ ° xx ^ dxdy ( X > V) 



Affine deformations of crystals on a cylinder 



u x (x,y) = a xx x + a xy y 
u y (x, y) = a yx x + a yy y. 



(84) 



These displacements will also result in a crystal, charac- 
terized by different unit vectors <?i and gi- The actual 
values of the matrix a will be determined by boundary 
conditions, such as counter-rotating twists applied to the 
ends of the cylinder. In general, the deformation matrix 
will not be symmetrical, although (neglecting a weak de- 
pendence of the lattice orientation relative to the cylinder 
axis) , the energy has to be independent of the antisym- 
(79) metric part \{a xy — a yx ). From the periodicity in the y 
direction we find that: 

u x (x, y + W) = u x (x, y) + b(mgi + ng 2 ), (85) 

where m and n are integers, leading to: 



W 



(mgi, x +ng 2 ,x)y Oty 



w 



{mg lty +ng 2 , v )- (86) 



Note that for small stresses corresponding to integers 
to, n which are not too large, we can replace the basis 
vectors g\ and g 2 by e\ and e 2 , up to corrections of order 
(b/W) 2 , which is a higher order effect. 

From Eq. ( 86 1 we conclude that the components of 



the deformation tensor a xy and a yy are quantized on 
a cylinder, while the components a xx and a yx are not 
quantized. Considering the strain matrix Uij = -^(diUj + 
djUi), we find that u xx = a xx is not quantized, since we 

a yy is quantized. 



have no restrictions on 



while 



L yy 



';!■:/ 



— 2 ( a xy + 



The off-diagonal component 
combination of a quantized object and a non-quantized 
object, and is thus non-quantized. The stress tensor is 
related to the strain tensor via the relation: 



= + \5ijiikk. 



(87) 



Each of the components of the stress tensor has a 
non-quantized piece. Thus, the Pcach-Koehler stress is 
fact non-quantized. Note, however, that particular lin- 
ear combinations of the stress tensor components can be 
quantized, similar to the case of superfluid films on cylin- 
ders pum]. 



Consider a lattice on a cylinder defined by minimal 
unit basis vectors e\ and e 2 . As discussed in section [J the 
lattice can be characterized by a pair of integers (M, N) 
such that: 



Mbei + Nbe 2 = Wy. 



(83) 
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